Parity and valley degeneracy in multilayer graphene 
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We study spatial symmetry in general yl_Byl-stacked multilayer graphene to illustrate how elec- 
tronic spectra at the two valleys are related in a magnetic field. We show that the lattice of 
multilayers with an even number of layers, as well as that of monolayer graphene, satisfy spatial 
inversion symmetry, which rigorously guarantees valley degeneracy in the absence of time-reversal 
symmetry. A multilayer with an odd number of layers (three or more) lacks inversion symmetry, 
but there is another transformation imposing an approximate valley degeneracy, which arises be- 
cause the low-energy Hamiltonian consists of separate monolayer like and bilayerlike parts. We show 
that an external electrostatic potential generally breaks valley degeneracy in a magnetic field, in a 
markedly different manner in odd and even multilayers. 

PACS numbers: 73.22.Pr 81.05.ue,73.43.Cd. 
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I. INTRODUCTION 

The fabrication of individual graphene flakes [l|, fol- 
lowed by an observation of the integer quantum Hall ef- 
fect in them [2,-4], triggered an explosion of interest in 
the electronic properties of graphene. It was fuelled, in 
part, by the realization that the low-energy band struc- 
ture of a graphene monolayer consists of two Dirac cones 
centered at inequivalent corners of the Brillouin zone, Kj^ 
and which are called valleys [Fig. [I](b)]. They sup- 
port chiral quasiparticles with opposite chirality in each 
valley, and a linear dispersion reminiscent of the quan- 
tum electrodynamics of massless fermions [Sj-tlj]. In the 
presence of time reversal symmetry, the energy spectrum 
is degenerate between the different valleys since the time 
reversal operation connects electronic states at to 
those at K^. In graphene, not only time reversal sym- 
metry, but parity, i.e., spatial inversion symmetry with 
respect to the center of a hexagon is also able to 

transform electronic states between valleys. In the pres- 
ence of a magnetic field, parity ensures degeneracy of the 
electronic spectra at different valleys. 

In this paper, we study spatial symmetry in general 
AB^-stacked (Bernal) multilayer graphene composed of 
N layers, to illustrate how the electronic spectra at 
the two valleys are related in a magnetic field. In 
multilayers with even N , including bilayers 0, [l^, [HI, 
the lattice obeys spatial inversion symmetry {x,y,z) — >■ 
(—x,—y,—z) similarly to monolayers, which swaps elec- 
tronic states between valleys ensuring valley degeneracy 
in the absence of time-reversal symmetry. The picture 
is different in multilayers with odd N, starting from 
trilayer graphene {N = 3), because their lattices do 
not satisfy spatial inversion symmetry Fig- 1(c). 

Here we consider an additional transformation which 
imposes an approximate valley degeneracy on the elec- 
tronic spectra of odd-A^ multilayers. It arises because 
the electronic Hamiltonian may be dec omp osed into sep- 
arate monolayer- and bilayerlike parts |13l - [l6| , with each 
part satisfying an approximate inversion symmetry cor- 
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FIG. 1: (a) Atomic structure of ABA-multilayer graphene. 
(b) A schematic of the Brillouin zone with two inequivalent 
valleys K± (c) Schematic of the lattice of bilayer graphene 
(left) and of Bernal-stacked trilayer graphene (right), with the 
inversion center for bilayer and the mirror plane for trilayer, 
respectively. 



responding to that of the real lattice of monolayer or bi- 
layer graphene, respectively. We predict a peculiar Lan- 
dau level spectra in bilayer and trilayer graphenes in the 
presence of interlayer potential asymmetry, with unusual 
structures of broken valley degeneracy that are markedly 
different from each other. 
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II. SYMMETRY ANALYSIS OF THE 
ELECTRONIC STRUCTURE 

We present here general symmetry arguments without 
referring to model-specific details, in order to provide in- 
formation about the valley degeneracy of Landau levels. 
We consider Bernal-stacked multilayer graphene with A^- 
layers in an external uniform magnetic field B with ar- 
bitrary direction, and an external potential U{r) (other 
than the lattice potential) where r — {x,y,z). Without 
resorting to any approximations, we can formally express 
the total Hamiltonian of the system as H [B, U{r)], and 
we take x and y-axes to be parallel to the graphene layers, 
and the z axis in the perpendicular direction. 

Even-N multilayers: In multilayers with an even num- 
ber of layers, the lattice is symmetric with respect to 
spatial inversion symmetry P [{x,y,z) — >■ (— a;, — ?/, — z)] 
as shown in Fig. [1] (c) , because the point group of the 
lattice Dsd [i,ll2 {{E,2C3,3C^,i,2S6,3ad}) can be re- 
garded as a direct product of group D3 ({£^, 2C3, 3C2}) 
with the inversion group Ci {{E,i}). In the absence of 
B, there is a symmetry with respect to time reversal T 
{t — > —t). The influence of P and T on the Hamiltonians 
of even-A^ multilayers is summarized as 

P: PH[B,U{r)]p-' = H[B,U{~r)] , (1) 
T: H*[B,U(r)] = H[-B,U{r)]. (2) 

The magnetic field B does not change sign under the 
operation P because it is an axial vector. Manes et al [9| 
used the combined PT operation to show that the Fermi 
points are stable with respect to the opening of a gap. 

In multilayer graphenes, electronic properties are well 
described by a k • p approximation in the vicinity of 
and points [I4I . We can consider the k • p Hamilto- 
nian (and thus eigenstates) at each valley separately, as 
long as the potential I7(r) is smooth in the a;?/-plane com- 
pared to the atomic scale, and the magnetic field is not 
too strong as to mix states around the two valleys. When 
B = 0, the time-reversal symmetry in Eq. ([2|) ensures 
valley degeneracy of the electronic spectrum, because the 
operation T swaps eigenstates at Kj^ and those of 
through flipping the Bloch factor as (e*^± '')* — e^^T "". 
Here K^(^ = ifc) are the wave numbers corresponding 
to points [141 • The operation of spatial inversion P 
also exchanges Kj^ and if_ because the point — K+ is 
equivalent to K_ in the Brillouin zone. Then the sym- 
metry of Eq. (HI) suggests that the eigenstates of at 
[B, U{y)] are related to those of at [B, U{—y)]. From 
this we immediately conclude that the Landau levels in 
even-A^ multilayers are degenerate in valleys as long as 
U{y) = U{—r), or the external potential has inversion 
symmetry with the same symmetry point as that of the 
lattice potential. 

Odd-N multilayers: In odd-A'^ multilayers with A" > 3, 
the point group i?3ft ({£;, 2C3, 3C^, a,,, 2^3, 3a„}) iHl 
can be regarded as a direct product of the group 
({i;,2C3,3C^}) with the reflection group Cs {{E^au)). 



With respect to the even- A" multilayers, spatial inversion 
is replaced by mirror reflection tr^ [(x, y, z) (x, y, — z)] 
as shown in Fig. [1] (c), which does not reverse the in- 
plane electronic momentum. The influence of T and ah 
on the Hamiltonians of odd-A^ multilayers is summarized 
as 



ah : ahH [B, U {x, y, z)] a,^ ^ 
T: H*[B,U{r)] 



H[B,U{x,y,^z)] ,(3) 
H[-B,U{r)]. (4) 



Time reversal symmetry, Eq. again ensures valley 
degeneracy of the electronic spectrum in the absence of 
an external magnetic field. The mirror refiection sym- 
metry would seem to play the role of parity. However, 
it does not transform between states at the two valleys, 
and is, therefore, unable to guarantee valley degeneracy 
in a magnetic field. It merely ensures that the spectrum 
at each valley is identical when the potential U (r) is in- 
verted with respect to z = 0. Actually, we can show 
that Landau levels in odd-layered multilayer graphenes 
are approximately valley-degenerate when U{r) = 0, by 
employing the effective mass Hamiltonian described in 
the following section. 

Although we focus on ABA graphene multilayers in 
this paper, we point out that ABC (rhombohedral) 
graphene multilayers, which have a different layer stack- 
ing 18, 19], have inversion symmetry irrespective of their 
layer number, and thus the valley degeneracy of Landau 
levels is always guaranteed. 



III. EFFECTIVE MASS HAMILTONIAN 

To investigate the Landau level structure in detail we 
adopt the effective-mass description of graphite in the 
Slonczewski-Weiss-McClure parameterization [2^. We 
consider AB-stacked A^-layer multilayer graphene in an 
external uniform magnetic field B perpendicular to the 
layer, and with an external electrostatic potential Uj at 
the j-th layer, which is uniform in the in-plane direction. 
In a basis with atomic components ipAi, '4'Bi, V'A2, V'-B2, 
''pASfJ^B^ • •ji the multilayer Hamiltonian at the val- 
ley mi is 



fHi V 

H2 V'' 
V V 

Hi T/t 



V 



(5) 



with 



V = 




{i: odd) 
(i: even) 



(6) 
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Here tt = (,Trx + ii^y, ^ = ±1 is the valley index and 
TT = —ihS/ + eA with the vector potential A, which 
gives the external magnetic field as B = V x A. Pa- 
rameter Ui is the potential at the i-th layer, which is 
constant within each layer, and S represents the energy 
difference between sites which have neighboring atoms 
right above or below them and those sites which do not, 
and thus it only exists for > 2. Parameter v is 
the band velocity of monolayer graphene, which is writ- 
ten as w = \/3ajo/2h. Other velocities are defined as 
vs = V^cL'-f3/2h and U4 = \/3a^4/2h. We neglect param- 
eters 72 and 75, which describe hopping between next- 
nearest neighboring layers. They actually break exact 
valley degeneracy as mentioned later. 

For even-iV multilayers with Uj = 0, the Hamiltonian 
Eq. ([5]) satisfies 



with (T2N being a 2N x 2N matrix. 



(7) 
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(8) 



which is nothing but the inversion symmetry discussed in 
the previous section. In odd-iV multilayers with Uj — 0, 
the Landau levels are degenerate in valleys due to a dif- 
ferent symmetry. We can show this by employing a uni- 
tary transformation which decomposes the Hamiltonian 
of A'^-layer graphene, Eq. ([5]), into subsystems equivalent 
to monolayer and bilayer graphenes. We construct the 
basis as [11 [H 



i=i 



where X = A,B, and 
2 



gniU) = - 



VN+1 
2 

VnTT 



(?) 

/TT ^ 

^cos(^-j) sm 



: sm 



cos 



2{N + 1) 

2(77TT 



-J 



(9) 

(10) 
(11) 



where j = 1, 2, • • • , iV is the layer index. The label m is 
the subsystem index which ranges as 



1,3,5,. 
0,2,4,. 



N : 

N ■ 



even 
odd 



(12) 



The superscript such as (A, odd) indicates that the wave- 
function has an amplitude only on \Aj) with odd j's. 

When the Hamiltonian ([SJ with Uj = is written in 
the basis Eq. the matrix is block-diagonalized in each 



m. The case of to = is special in that gm{j) is identi- 
cally zero, so that only two bases {\(f>Q^'°'^'^'') , \4'(^'°'^'^^)} 
survive in Eq. Q . The submatrix is written for this two 
component basis as 







(13) 



which is equivalent to the Hamiltonian of mono- 
layer graphene except for the diagonal terms con- 
taining 5. For m 7^ 0, the submatrix for 



(A, odd) 



,(B,odd)i 



H 



(0 



(A,cvcn)\ I ,(B,cvcn)\T , 

^ ''^ )| becomes 

(14) 



VTT S A71 

Xv^TT A71 6 

y AwaTT^ —\v4Tr vtt 



-Aw47rt 

rt 



VTT' 





where A = Am is defined by 



2(iV+ 1)' 



(15) 



Eq. (|T4|) is identical to the Hamiltonian of bilayer 
graphene, except that interlayer-coupling parameters 71, 
73 and 74 are multiplied by the factor A. 

Since the decomposed Hamiltonian matrices are analo- 
gous to those of monolayer or bilayer graphene, they obey 
the corresponding inversion symmetry. The bilayer-type 
submatrix, Eq. ((Tl)) . obeys 



(16) 



which guarantees valley degeneracy of Landau levels. 
In the original basis, this "effective" inversion process 
Eq. (jl6p exchanges the wave amplitudes on odd-jth lay- 
ers and those on even-jth ones. In even-iV multilayers, 
the operation Eq. ((T6| for each eigenstate becomes equiv- 
alent to the original inversion symmetry Eq. ([7]) except 
for a phase factor. For the monolayer- type submatrix Eq. 
([T^ (exists in odd-A^) satisfies 



^2 ^m=0 ^2 ~ n„i=' 



0' 



(17) 



but it is only approximate since S breaks this symmetry, 
unlike in the bilayer-type symmetry, Eq. (jl6p . The oper- 
ation of Eq. (I17[) exchanges the wave amplitudes within 
odd-j layers. 

The extra parameters 72 and 75 neglected here gener- 
ally mix the states between different to's, and also appear 
in diagonal elements within each m [l^j. In odd- TV mul- 
tilayers, they lift the valley degeneracy by breaking the 
effective inversion symmetry Eq. (jl6p through additional 
matrix elements. It should be noted that, in even-iV mul- 
tilayers, valley degeneracy is never infiuenced by any ex- 
tra lattice parameter, because it is protected by inversion 
symmetry inherent in the lattice. 



4 




FIG. 2: The low-energy Landau level spectrum of bilayer 
graphene with rfl/71 = 0.3 at each valley, plotted as a func- 
tion of interlayer asymmetry A. 



IV. LANDAU LEVEL SPECTRA OF BILAYER 
AND TRILAYER GRAPHENE 

Here, we illustrate the valley degeneracy of the elec- 
tronic spectra in even- and odd-A^ multilayers through an 
analytical description of the Landau level spectra of bi- 
layer and trilayer graphene. We use the magnetic length 
Xb = \/ h/{eB) and an energy scale Tb — \/2hv^eB — 
\/2hv/\B related to the inverse of the magnetic length. 
The Landau level spectrum, in a magnetic field of mag- 
nitude B perpendicular to the graphene sheet, may be 
found using the relation (tt, tt^ ) = {y/2h/XB){a\a) for 
K+, and (7r,7rt) = {V2h/ XB){a, a'<) for where 
a' and a are raising and lowering operators, respectively, 
which operate on the Landau-level wave function ipn as 
a(pn = y/nifn-i and a^n = Vn + lifn+i- 

Bilayer graphene: the Hamiltonian of bilayer graphene 
[Tsl . [24 . [2H is given by the multilayer Hamiltonian, 
Eq. ([S]), with N = 2. Here, we consider a simple form 
of the Hamiltonian that only contains non-zero param- 
eters 7o, describing nearest-neighbor intralayer hopping, 
71, describing the dominant interlayer coupling, and the 
layer potential defined as {1/1,1/2) = A(l,— 1). In the 
absence of interlayer asymmetry, A = 0, the spectrum 
consists of fourfold (valley and spin) degenerate Lan- 
dau levels [l0|, except for the level at zero energy which 
is eightfold degenerate. The valley degeneracy is guar- 
anteed by spatial inversion symmetry Eq. ([1}. Finite 
A breaks spatial inversion symmetry and it splits the 



valley degeneracy of the levels [10|, I26l - l30l | , as indicated 
in the numerically-calculated spectrum with a magnetic 
field of rB/71 = 0.3, plotted in Fig. M [U. There we 
can see that the energy levels of two valleys are related 
as eif^(A) = A), due to the inversion symme- 

try arguments presented in Sec. HIl The valley splitting 
Ek^ (A) — Ek^ (A) is thus an odd function in A, and gen- 
erally begins with a term linear in A. 

In the limit {E,^/n^B,^} ^ 71, the Hamiltonian is 
approximately described by an effective Hamiltonian (Toj 
operating in the space of two-component wave functions 

IpAl, i>B2 



n 



(off) 
AB 







+A 



- % 
71 Vtt^ 





-1 + 2v'^7nTy-ff 



.(18) 



The energy levels are obtained by assuming the wave 
function '9n,K^ = {citpn+^,C2ipn-0, whcrc (pm<a is re- 
garded as 0. They are given by 



£n>i — 



£0 
£-1 



.API 

,a(i-^ 



A2 



1 - 



2r| 

7? 



(19) 

(20) 
(21) 



The splitting of the valley degeneracy in n > 1 (by terms 
containing ^) is linear in the asymmetry potential A, and 
also in the magnetic field B. 

Trilayer graphene: The trilayer graphene (l2| - [T5l [2ll . 
|32| - [36| Hamiltonian is given by the multilayer Hamil- 
tonian, Eq. ([5]), with = 3. We again retain the 
parameters 70 and 71, and set the layer potential as 
(1/1,1/2,1/3) = Ai(l,0, — 1), to focus our attention on in- 
terlayer asymmetry Ai — [Ui — U3) /2. As described 
in Section IIIH we perform a unitary transformation to 
decompose the Hamiltonian into a monolayerlike part, 
Eq. (IT51) . and a bilayerlike part, Eq. (ITil) (such a de- 
composed trilayer Hamiltonian is written explicitly in 
Ref. [131). 

In the absence of interlayer asymmetry, Ai =0, the 
spectrum consists of superimposed monolayerlike and bi- 
layerlike spectra [13| , with fourfold (valley and spin) de- 
generate Landau levels, except for the level at zero energy 
which is twelvefold degenerate. The valley degeneracy is 
guaranteed by the effective spatial inversion symmetry 
of the separate, monolayerlike and bilayerlike, parts of 
the decomposed Hamiltonian. Finite Ai breaks the ef- 
fective spatial inversion symmetry and it splits the valley 
degeneracy of the levels as indicated in the numerically- 
calculated spectrum at rB/71 = 0.3, plotted in Fig. [3] 
[sij. Unlike in bilayer, the energy spectrum at each val- 
ley is an even function of Ai, due to the reflection sym- 
metry argued in Sec. [Hi The valley splitting is therefore 
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For index n = 0, there are three energy levels, 



FIG. 3: The low-energy Landau level spectrum of ABA- 
stacked trilayer graphene with Vb/'^i = 0.3 at each valley, 
plotted as a function of interlayer asymmetry Ai. 



quadratic in Ai, except for zero energy where the differ- 
ent levels are degenerate at each single valley, and the 
energy level (and thus splitting) is allowed to be linear 
in Ai. 

In the limit {e, v^Fs, Ai} ^ 71, it is possible to ob- 
tain a simplified description of four electronic bands near 
zero energy, by eliminating atomic components associ- 
ated with bands split away from zero by energy ±'\/27i, 
to obtain a four-component effective Hamiltonian in ba- 
sis [ll}Al - ilA3\/V2, [ipBl - l\jB-i\l\/^, {^A\ + ■(/'A3]/V2, 
V'B2, 



n 



(eff) 
ABA 



( 

VK 

Ai 

V 





Ai 














Aii;7r 




\/27i 


V27I 



\ 

\/27l 
\/27l 

/ 



(22) 



The energy levels are obtained by assuming the wave 
function ^n,K^ = (ci<^„+4, C2(y5„, C3<^„+^, C4<y5„-4). For 
levels with index n > 1 at valley 7^+, the Hamiltonian 
Eq. (j22p yields four energy levels for each n, given by 



-(") 

'n>l,K+ 



= ± 



V^Fb [Aj-(n + l)F|] 



(23) 
(24) 



Jl) 



±a/A?+F|, 



= 0, 



and, for n = — 1, two energy levels. 



-(«) 



±Ai 



(25) 
(26) 

(27) 



Here, the superscript (u) means that the level is related 
to an "upper" energy band which is split away from zero 
energy by Ai (i.e. the Landau level energy approaches 
±Ai as Tb — 0), whereas superscript (l) means that 
the level is related to a "lower" energy band near zero 
energy (i.e. the Landau level energy approaches zero as 
Tb^O). 

The Landau level spectrum at the second valley, 
differs as compared to that at . For levels with index 
n > 1 , the Hamiltonian Eq. ((22)) again yields four energy 
levels for each n given by 



'n>l,K' 



-n>l.K- 



^ ±JAl + nrl 



(28) 



7271 y/A'i + nTl ■ 
For index n ^ 0, there are two energy levels at 

AiFb 



-(0 

'-0.K- 



±- 



V271 

and, for n — —1, one energy level. 



.(0 

'--l.K- 



0. 



(30) 



(31) 



In the limit Ai — 0, the levels are degenerate in pairs of 

= (" > 0)' and = {n > 1), 

which form monolayerlike and bilayerlike spectra, respec- 
tively. In the high-field limit, |Ai| ^ F^, the valley 
splitting of those levels is proportional to Af /71 but in- 
dependent of B, unlike bilayer graphene. 

The other levels give degenerate zero-energy states at 

Ai = 0. For K+{K^), one of the levels £^"^+(4!^-) 
corresponds to the zero-energy level of the monolay- 
erlike spectrum, whereas ^^_|.(e!^i /f_) and one of 



e^^l j^j^{eQ'j^_) correspond to two zero-energy levels of 
the bilayer-like spectrum. This yields an overall twelve- 
fold degeneracy with spin degeneracy included (as op- 
posed to the eightfold degeneracy in bilayers and fourfold 
degeneracy in monolayers). In the presence of finite inter- 
layer asymmetry, two of the otherwise-zero levels at each 
valley are hybridized and split away from zero [Eq. (|27p at 
and Ea. (l5D|) at i^-], whereas one level at each valley 
remains at zero [Eg.pG]) at and Eg. ([31]) at K-\. Zero 
energy states are formed primarily by different atomic or- 
bitals at different valleys, and, as in bilayers fiol, ITsl. [37| . 
the orbitals in trilayers are on different layers, so that 



6 



those levels exhibit quite different dependences on the 
intcrlayer asymmetry between and K_ . 

In realistic experimental situations, inter layer poten- 
tial asymmetry can be produced by gate-induced electric 
fields. There the potentials Ui, which are taken as ex- 
ternal parameters in the present work, should be deter- 
mined self-consistently including the screening effect of 
graphene electrons ^24, 32]. The self-consistent Landau 
level structure and quantum Hall effect as a function of 
gate voltage would be an important future study [s^ . 

It has been pointed out [38| that a number of level 
crossings will occur in bilayer graphene at finite asym- 
metry. Fig. [2] We note that the Landau level spectrum 



of trilayer graphene. Fig. [3l has a very rich pattern of 
level crossings at a range of different values of energy 
and asymmetry. The presence of additional terms in the 
Hamiltonian will tend to produce some anti-crossings in 
the spectrum. The precise position and nature of level 
crossings in the spectra of multilayer graphenes will be 
the subject of future investigation. 
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